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We study numerically the behavior of a single quantized vortex in a rotating semi-infinite cylinder. We study 
in particular the spiraling motion of a vortex in a cylinder that is parallel to the rotation axis. We determine 
the asymptotic form of the vortex and its axial and azimuthal propagation velocities under a wide range of 
parameters. We also study the stability of the vortex line and the effect of tilting the cylinder from the rotation 
axis. 



I. INTRODUCTION 



Turbulence is a much studied, but not yet fully understood, 
phenomenon in both classical and quantum fluids. It is hoped 
that research in quantum turbulence will also shed light on 
classical turbulence, besides being an interesting field of re- 
search in itself 1 ' 2 . To study quantum turbulence one needs to 
study the motion of quantized vortices. This has been done 
using various experimental and theoretical approaches 1 ' 3-5 . In 
theoretical approaches the emphasis has been shifting to nu- 
merical simulations. Many of these include a large number 
of vortices and require a considerable amount of computing 
power, especially when calculations are performed using the 
full Biot-Savart equations 6 . However, there is still a need to 
better understand the motion of a single vortex. Since the mo- 
tion of a curved vortex line is somewhat counter-intuitive, and 
solving the equations analytically is difficult and prone to er- 
rors, it is convenient to use numerical simulations. For single 
vortices, it also becomes a realistic possibility to scan a large 
volume of parameter space, i.e. various combinations of pres- 
sure p, temperature T, rotational velocity ft, vessel size and 
shape, and initial vortex configurations. 

An existing computer software, previously used mainly for 
studying the large scale behavior of many vortices 3 ' 7 ' 8 , is ap- 
plied to study the motion of a single vortex line. More specif- 
ically, we study the motion of a single superfluid vortex fil- 
ament in a rotating semi-infinite cylinder, as illustrated in 
Figs. 1 and 2. Earlier studies by other researchers have con- 
sidered a single vortex in a cylinder attached to a wire 9 ' 10 . 

We assume the normal component of the velocity to be in 
solid-body rotation. Our main point of interest is the asymp- 
totic velocity of the vortex end touching the side wall of the 
cylinder. We calculate the evolution of the vortex from the 
initial state using the vortex filament model of the two-fluid 
paradigm with Biot-Savart formalism, until the asymptotic 
situation is reached. Although we have used 3 He-B -specific 
parameters in our calculations, the methods and the results 
can be generalized to vortices in other superfluids and Bose- 
Einstein condensates in many cases. 

One surprising effect discovered in our calculations is the 
stability of the vortex, even in the low temperature limit. This 
contradicts the expectations based on some earlier simulations 
and experimental results in 3 He-B 11-14 . We assume that the 
discrepancy with experiments is caused by some surface ef- 
fects not accounted for in our model (such as pinning), or due 




FIG. 1 : The evolution of a vortex in a rotating coordinate system. The 
total time span 2.79/aQ comprises 61 snapshots. Other parameters 
are Re a = 3.63 (~ 0.4T C in 3 He-B) and QR 2 / 'k = 13.6. 



to some uncontrolled heat leaks or superflows in the experi- 
ments that may destabilize the vortices. 



II. THE MODEL 

We study a superfluid in a half-infinite cylinder of radius 
R that rotates at constant angular velocity tt. The axis of 
the cylinder may be tilted by an angle 6 relative to the rota- 
tion axis. We use cylindrical coordinates (r, 0, z) fixed to the 
cylinder. Our study is based on the two-fluid model, where 
the normal and superfluid components have velocities v n and 
v s , respectively. The normal component is assumed to be in 
rigid body rotation, v n = ft x r. This is a good assumption in 
3 He-B, because mutual friction is weak in comparison to the 
viscosity of the normal fluid 15 . 
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b) 270 

FIG. 2: The evolution of the vortex shown in Fig. 1 presented in 
cylindrical coordinates fixed to the rotating frame, a) (r, z) plot with 
37 snapshots in time span 4.19/aQ. b) Polar plot with 31 snapshots 
in time span 1.40/aQ. 

The superfluid velocity v s is determined by vortex lines 6 . 
The location of the vortex lines are given by s(£, t), which is 
parametrized by the vortex length £ and time t. Partial deriva- 



tives with respect to £ are denoted by primes. Then, the unit 
tangent of the vortex core line is the first derivative of s with 
respect to £. This is denoted by s' (or s' to emphasize that it is 
a unit vector). The superfluid velocity v s is calculated using 
the (full) Biot-Savart formalism: 

v s (r,t) = v w + v b , (1) 

with 

_ k_ f (si - r) x rfsi 

where C denotes all the vortex lines and k is the circulation of 
the superfluid velocity around a vortex line; all the positions 
and velocities are in laboratory coordinates. The first term, 
v w , in (1) needs to be numerically calculated from the integral 

(2) . It has a divergence on the vortex core, which has to be cut 
off at the vortex core radius a, see Refs. 6 and 16 for details. 
The second term, Vb, is the boundary field, or image velocity 
field, needed to prevent flow across the vessel boundary. The 
boundary field v b can be calculated from the Laplace equa- 
tion, or by the method of image vortices. To calculate Vb we 
use the method of approximate image vortices as in Ref. 16. 

The derivative of s(£, t) with respect to time defines the vor- 
tex line velocity vl except that the component of vl parallel 
to the line is not defined. The equation of motion is commonly 
written in laboratory coordinates as 6 

v L = v s + as' x (v n - v s ) - a's' x [§' x (v n - v s )]. (3) 

It contains the mutual friction parameters a and a' . These 
parameters for 3 He depend (dramatically) on temperature and 
(moderately) on pressure 17,18 . Whenever not explicitly speci- 
fying their values, we use their temperature dependent values 
for 3 He-B at the pressure of 29.34 bar. 

Using the fact that the normal component is in rigid body 
motion (v n = ft x r), we can rewrite the equation of motion 

(3) in rotating coordinates (where v n = 0) as: 

v L = a(Re a v s - s' x v s ). (4) 

Here we have defined Re a := (1 — a')/ a, which has some 
similarity with the Reynolds number defined in normal fluids. 
In 3 He-B the quantity Re a (T) decreases monotonically from 
oo to 0, when the temperature increases form to T c . From 
now on all velocities will be in rotating coordinates, unless 
explicitly stated otherwise. 

We study in particular the velocity of the vortex end in the 
long-time limit (t — » oo). For that we define vl z and vl(/> by 
writing in the rotating frame 

v L v Lz e z +tf£,0e0. (5) 

Note that typically VL<f> and ft have opposite signs: for a grow- 
ing vortex, as in our case, vlq < 0, if Vt > 0. 

It is also possible to define the climbing angle j3 of the vor- 
tex in the rotating frame by tan f3 = vlz/vl(/>- The climbing 
angle in the laboratory frame satisfies tan/3 = VL Z /(vLcf) + 
SIR). 



III. SCALING PROPERTIES 

The parameter space of the system (without considering the 
parameters describing the initial vortex configuration) con- 
sists of mutual friction parameters a and a! , the vortex core 
radius a, the circulation quantum k, the cylinder radius R, the 
angular velocity Q, and for a tilted cylinder the tilting angle 

0. The exploration of this parameter space is simplified by 
dimensional analysis. Let us, for example, study the asymp- 
totic axial velocity vl z (5) of the vortex end in the absence of 
tilting, = 0. According to dimensional analysis, the dimen- 
sionless velocity vl z /VlR can only depend on dimensionless 
quantities appearing in the problem, which in this case are 
R/a, VtR 2 1 'k, a and a'. Thus the parameter space is four- 
dimensional. Using properties specific to the present prob- 
lem, we show in the following that this space can further be 
reduced to three dimensions exactly and to two dimensions 
approximately. 

Let us consider the evolution of a vortex state from a fixed 
initial configuration. We see from Eq. (4) that varying a but 
keeping Re^ constant has the effect of changing the time scale 
only. This means that the dimensionless time is aQt and the 
dimensionless line velocity is v^/atiR, but otherwise there is 
no dependence on a. (Note that the dimensionless superfluid 
velocity v s /QR does not have a.) Since the asymptotic ve- 
locity is obtained in the limit t —> oo, the dependence on a fit 
drops out. Since the system is dissipative (for a ^ 0), the 
same final state can be obtained by a variety of initial states, 

1. e. the details of the initial state are not important either. Thus 
we conclude that the asymptotic velocity has the form 

/OR 2 R l-a f \ 

v Lz aQRf , -, (6) 

\ k a a J 

with some function /. The same form applies to the azimuthal 
velocity vlq- Note that this is valid only in the rotating frame; 
the laboratory frame azimuthal velocity vl^ + QR does not 
have the form (6). Because the parameter VtR 2 / 'k appears 
repeatedly, we define a unitless rotational velocity: 



We also remind that the last parameter in (6) has the short- 
hand Re a = (1 — a!) /a. 

An approximate simplification is offered by the fact that 
vl has only a weak logarithmic dependence on R/a, see 
e.g., Refs. 19 and 20. Hence, different values of R/a have 
only a minor effect on the the results; this is demonstrated 
in Fig. 3. Because of this weak dependence, we have fixed 
R/a = 1.77 x 10 5 in all other calculations in this paper ex- 
cept the ones presented in Fig. 7. 

IV. ASYMPTOTIC VORTEX VELOCITY IN A 
NON-TILTED CYLINDER 

The vortex evolution is solved by numerically integrating 
Eq. (4) forward in time in the rotating frame. We use a quar- 
ter vortex ring as our initial condition. The radius of the ring 
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FIG. 3: The effect of the parameter R/a. The solid lines have R/a 
equal to 1, 2 and 4 times 1.77 x 10 5 while other parameters have 
constant values (aQt = 0.75, Re a = 3.63 and v = 13.6). These 
can be contrasted to two cases, which differ by ±10% in the value 
of Q from the middle case (the dashed line has aQt = 0.82 and 
v = 15.0 and the dash-dotted line has aQt = 0.67 and v = 12.2). 
The vortices are shown in the rotating coordinate system. 

is half of the cylinder radius. The ring is situated so that the 
vortex starts from the bottom of the cylinder and ends at the 
cylinder wall. Both ends of the vortex are allowed to move 
freely. Any vortex end touching the wall is always normal to 
it. In principle, a vortex line may also form a closed loop, 
without touching the boundary at all. A vortex loop may be 
created when a vortex line reconnects with itself. If a vor- 
tex loop comes very close to the wall, it reconnects with it, 
i.e. snaps open, with the two vortex ends connecting to the 
wall orthogonally. The use of a quarter vortex ring as the ini- 
tial condition is not necessary, but is convenient in our case. 
What matters is whether the required asymptotic evolution (if 
it exists) is reached from the initial condition or not, and a 
quarter ring is suitable in a large parameter range. However, 
if the temperature is low and the shape of the vortex is far 
from the asymptotic one, the vortex may oscillate a long time 
before reaching a stable form. For Re^ > 100, even the quar- 
ter ring vortex leads to long-lasting oscillations in the vortex 
shape. 

Our goal is to study the velocity of the end point of the 
vortex vl (5). Before going to the general case, we consider 
two limiting cases where analytic solutions are found. 

In the limit v ^> 1, the cylinder wall at r = R can be 
thought of as a plane and we suppose that the vortex line ap- 
proaches the wall as a straight line normal to it 27 . We can now 
put ft = Qe z , r = Re r , a point on the large cylinder wall 
(or plane) and s' = e r . This gives v s = —flRe^ in rotating 
coordinates, and inserting it into (4) we get: 

v L = SlR[ae z + (a' -1)$,/,]. (8) 

Hence vl z = OiQR, vl^ = (a' — l)QR, and tan (3 = —Re" 1 . 
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FIG. 4: The axial velocity of the vortex end vlz, plotted as a func- 
tion of v = QR 2 /k at Re a values 78.6, 20.6, 7.65, 3.63, 2.18, 1.46, 
1.03, 0.74, 0.53, 0.36, 0.24, 0.15, and 0.083, corresponding to tem- 
peratures T/T c = 0.25, 0.30, 0.35, 0.85 in 3 He-B, respectively, 
a) Plot in linear scale, b) Logarithmic plot. 



FIG. 5: The azimuthal velocity of the vortex end vl<i>, plotted as a 
function of v at the same values of the parameter Re a as used in 
Fig. 4a and Fig. 4b, again corresponding to temperatures T/T c = 
0.25, 0.30, 0.35, 0.85 in 3 He-B. a) Plot in linear scale, b) Loga- 
rithmic plot. 



Another exact result corresponds to vanishing velocity, 
vl = 0. This is obtained in the case that the one- vortex 
and no- vortex states are in equilibrium. This means that there 
is no force to drive the vortex and therefore the vortex end 
must be stationary in the rotating frame. The critical rota- 
tional velocity corresponding to this will be derived below, 
and is given in Eq. (16). In the dimensionless form it corre- 
sponds to v c = Q c R 2 /k = (l/27r) ln(R/a), which in our 
case (R/a = 1.77 x 10 5 ) has the value 1.92. 

The numerically calculated vl z is depicted in Fig. 4a. It 
shows vlz as function of v at different temperatures. The low- 
est point corresponds to v c , where vl z vanishes. For large v 
the result tends to the analytic result (8). These two limits 
are obeyed at all temperatures, but a temperature dependence 
remains in between. For quantitative purposes it seems rea- 
sonable to write 



VLz 



aSlR(l - e z ). 



(9) 



where e z — >• for v — >• oo. The value of e z defined hereby is 
shown as a function of v as a log-log plot in Fig. 4b. 

Figure 4b shows at low values of Re^ (or, equivalently, at 
high values of T/T c ) an almost perfect linear relationship be- 



tween In e z and In v. In order to fit this, we keep the first 
(analytically obtained) point (y — v c , e z = 0) fixed: 



lne 2 



k z In - 



(10) 



With Re a = 0.083 (~ 0.85T C in 3 He-B), fitting in the range 
v e [i/ c , 871] gives k z ~ -0.741. Thus, 



VLz 



attR 



1 



(V_c_ 

V v 



0.741 



(ID 



We now proceed to analyze the azimuthal velocity v^. The 
numerically calculated vl^ is shown in Fig. 5a. Similarly as 
for VL Z /m order to resolve more details we define by 



^ = (c/-l)OR(l-e ). 



(12) 



The log-log plot of vs. v is shown in Fig. 5b. 

We see that there is an approximately linear dependence of 
In on In v. 



In e d 



krh In - 



(13) 
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FIG. 6: Asymptotic vortex forms plotted for combinations of low 
and high temperatures (0.25T C ; 0.95T C ) and rotational velocities {y c \ 
452.9z^ c ): a) r against z, b) <j) against z, c) (j) against r. The vortex 
shape at v — v c is independent of temperature, so only the solution at 
0.25T C is plotted. Another point to make is the negligible variation 
of 4> in all the four cases, except for the combination of high rota- 
tional velocity and low temperature (452.9z^ c , 0.25T C ), dash-dotted 
line. These solutions can be compared to the approximative (LIA) 
solution found in Eq. (13), Ref. 21, dotted line. The temperatures 
are approximative values for 3 He-B. The Re a parameters used in 
the calculations are: Re a = 78.6 for 0.25T C and Re a = 0.016 for 
0.95T C . 
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FIG. 7: Comparison of numerically calculated asymptotic velocities 
between a single vortex end and a vortex front. The lines are for a 
single vortex, giving the axial (solid line) and azimuthal (dash-dotted 
line) components. The data points are for the vortex front, giving the 
axial (circles) and the azimuthal velocity (crosses). All azimuthal 
velocities are in the rotating coordinate system. The front velocity is 
from numerical simulations, originally presented in Ref. 4, p. 3219. 
The parameters are v = 34.0 and R/a = 0.88 x 10 5 . 

At high temperatures the axial velocity (Fig. 4b) closely 
follows the formula in Eq. (11), whereas the azimuthal log- 
log plot, Fig. 5b, looks linear only when In v > 3. Indeed, 
fitting in the range In v G [3.30, 6.77] at Re a = 0.083 gives 
fc -0.719 with v ~ 1.107. 

At low temperatures, the opposite seems to be the case; 
the axial velocity is nonlinear in the log-log plot (Fig. 4b), 
while the azimuthal velocity (Fig. 5b) seems to be linear at 
In v < 4. Keeping the first point fixed and fitting in the range 
In ^ £ [In v c , 4.00] at Re a = 78.6 (0.25T C in 3 He-B) gives 
kff, = —0.98. (Fitting without fixing the first point gives 
k(f) = —0.99 and vo ~ 1.97.) At even lower temperatures 
the vortices converge very slowly to the asymptotic form, so 
it is difficult to say what the behavior is like at T = 0. 

The shape of the top part of a rotating vortex is depicted 
in Fig. 6. The asymptotic form depends somewhat on the pa- 
rameters. We see that the shape at the critical value v = v c is 
independent on temperature. This is natural since the vortex 
is at rest in the rotating frame and thus the dissipation does 
not play any role. The azimuthal angle of the vortices seems 
constant both at v — v c and at high temperatures. We also 
compare the shape with the one obtained using a simple model 
having an analytic solution 21 . 

We can take the expression (1 1) as a basis for further analy- 
sis, including e.g., the analysis of the axial velocity of a mov- 
ing vortex front. The comparison to the vortex front is made 
in Fig. 7. The velocity of a moving front is affected by the 
twisting of vortices behind the front. An analytic formula has 
been applied to analyze the effect of this twisting on the front 
velocity 3 ' 22 . What is studied in this paper could be termed 
as the single vortex contribution to the vortex front velocity; 
twisting of vortices adds another contribution, surface friction 
yet another, and phenomena associated with turbulence (re- 
connections, Kelvin cascade, vortex tangle diffusion, bound 
fermion states at vortex core) yet another. 
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V. STABILITY OF VORTICES IN NON-TILTED 
CYLINDER 

A. Our definition of stability 

We now consider the stability of vortices in a non-tilted 
cylinder. From the standpoint of this paper, stability means 
that the vortex configuration has an asymptotic state with a 
single vortex, moving upwards and rotating with a constant 
angular velocity without changing its shape. In order to reach 
this kind of asymptotic vortex configuration, the angular ve- 
locity ft must be above some critical value Q c . Otherwise, the 
vortex will shrink, and eventually disappear. 



P4 

a 



0.04 



0.03 



0.02 



0.01 



-0.01 




B. Critical angular velocity 

The critical angular velocity ft c depends on the cylinder ra- 
dius R. The necessity of the critical value can be heuristically 
justified by recalling that vortex rings need a sufficiently large 
counterflow being pushed through them in order to grow, and 
furthermore, realizing that the same is true for vortices that 
end on walls: counterflow needs to be pushed through the area 
delimited by the vortex line and the walls. The counterflow in 
this case is provided by the rotation of the vessel. Without a 
numerical solution, this idea only provides a rough approxi- 
mation of the value of Q c . 

An exact value for Q c , however, can be obtained by study- 
ing the free energy 23 



F — £?kin — QL, 



(14) 



where E^ n is the kinetic energy and L is the angular momen- 
tum of the fluid. In order for a one- vortex state to be stable, F 
should be smaller for one axial vortex than for the stationary 
superfluid. This gives the condition 



— In— - -Q(i? 2 - a 2 ) < 0. 



(15) 



We neglect the small term ex a 2 . Thus for the one-vortex state 
to be stable, ft has to be larger than 



Q r = 



k R 
m- 



2ttR 2 a ' 



(16) 



In dimensionless form this can be written v c = Q C R 2 / k = 
(l/27r)\n(R/a). 



C. Other causes for instability 

Our numerical calculations indicate that once we have a 
vortex, whose end moves upwards, it is found to be stable at 
all temperatures, even at the zero temperature limit. Instability 
is found only if the angular velocity is below the critical value 
Q c , or if the cylinder is tilted. This is somewhat surprising, 
since previous simulations and experimental studies in 3 He- 
B have indicated that below a certain temperature, roughly 



FIG. 8: Value of counterflow projected on the vortex tangent, non- 
tilted cylinder. There is a clear peak near the top of the vortex, then 
a rapid drop to zero. The critical value of projected counterflow ve- 
locity for Kelvin waves with wavelength R is v c ,og = 0.1634 QR, 
which is clearly not reached. The parameters are the same as for the 
last vortex state in Fig. 2a. 



0.5T C (or Re a > 1.5), single vortices become unstable 2 ' 11 ' 14 . 
We suggest that this kind of single vortex instability is due to 
some surface effects, or heat leaks causing extra counterflow. 
In our model, it is assumed that the surface of the cylinder is 
smooth and that there is no surface-induced extra friction at 
the vortex end. 

Another source of instability is related to the creation and 
growth of Kelvin waves. The so-called Ostermeier-Glaberson 
instability 24 ' 25 appears when the counterflow v c f := v n — v s 
along the vortex line is sufficiently strong, and causes straight 
vortex lines to become unstable towards the appearance of 
Kelvin waves. For a single axially oriented straight vortex, 
the critical counterflow velocity above which the amplitude of 
Kelvin waves with a wave vector k are able to grow is: 



fl + r]k 2 t k , 1 

^,og = ^— , where, = - In-. 



(17) 



In our case of a non-tilted cylinder with a vortex attached 
to the cylinder wall this effect may exist in some situations. 
From Fig. 2b one observes that the vortex has a small com- 
ponent along the azimuthal direction. This results in a non- 
zero counterflow along the vortex and could in principle cause 
Ostermeier-Glaberson instability. However, the numerics in- 
dicate that the Kelvin waves that could grow have a wave- 
length much longer than the cylinder radius (and therefore 
larger than the region of this finite counterflow). Hence, the 
Kelvin waves cannot grow and Ostermeier-Glaberson insta- 
bility does not play any role. There is in principle room for 
Kelvin waves to grow along the lower part of the vortex as 
well, but since the Kelvin waves that are created in the un- 
tilted case move upwards along the vortex line, and since the 
counterflow along the vortex line has its maximum in the top 
part of the vortex and falls sharply below it (Fig. 8), there is 
no chance for Kelvin waves to grow in the lower parts of the 
vortex either. However, the situation in a tilted rotating cylin- 
der is different. With enough tilting, the boundary superfluid 



FIG. 9: Kelvin waves created in a rotating tilted cylinder. The ap- 
pearance of Kelvin waves, (which can be clearly seen in the lower 
left quadrant of the polar plot) is due to Ostermeier-Glaberson insta- 
bility caused by the background velocity field Vb. This background 
velocity field is nonzero in a rotating tilted cylinder and in this case 
causes a sufficiently strong counterflow along the vortex line to cre- 
ate Kelvin waves. The parameters are aQt = 0.56, Re a = 20.6 
(~ 0.3T C in 3 He-B), QR 2 /k = 13.6 and = 60°. 



velocity Vb creates a sufficiently strong counterflow, which re- 
sults in the appearance of Kelvin waves. These Kelvin waves 
then produce a change from a laminar to a more turbulent fluid 
motion, possibly resulting in a vortex tangle, see Fig. 9. Yet 
another phenomenon is the Crow instability 26 which applies to 
vortices that are close to each other, or close to the wall (which 
can be interpreted as close to the image vortex). This instabil- 
ity is due to the fact that vortices of opposite direction attract 
each other. Thus, the parts of vortices that are slightly closer 
to each other start moving faster towards each other, deform- 
ing the vortices even more, and eventually resulting in a large 
number of reconnections and new vortices. In our case the 
Crow instability does not work in the untilted cylinder since 
the rotation provides a stabilizing effect. 

These various causes of instability notwithstanding, we 
have some doubts about the equations for vortices in a rotating 
container used in Refs. 11-13. According to our calculations 
the creation of turbulent states is more difficult than the results 
in those papers would suggest. 



D. Effects of initial vortex state 

While it is virtually impossible to carry out an exhaus- 
tive study for this case, since the possible space of initial 




FIG. 10: Vortex evolution for different tilting angles. Twenty-one 
snapshots are shown for each tilting with angles = 10°, 20°, 30°, 
40°, 50°, and 60°. With increasing tilting some qualitative differ- 
ences appear. The starting point of the vortex may climb on the side 
wall. The end point may start moving axially. All plots are in the 
rotating coordinate system. The total time is 2.33/aQ, Re a = 3.63 
and v — 13.6. 



conditions has an infinite-dimensional parameter space, some 
things can still be said. 

The initial vortex state used here was a quarter ring attached 
at the bottom and at the side wall. This was sufficient to 
guarantee the required asymptotic state. Before the vortex as- 
sumed a behavior close the asymptotic one, there were some 
transient effects. The velocity of vl z oscillated before reach- 
ing the asymptotic value. This oscillation was stronger at low 
T/T c . Furthermore, if the initial radius was large, this oscil- 
lation was usually stronger. This can be explained as follows: 
At small scales, the fluid can always be considered as a bulk 
fluid having uniform normal and superfluid velocities. For 
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these uniform velocities, the non-oscillating shape for a closed 
vortex is a perfectly circular ring. Also, at small distances, the 
curvature of the walls can be neglected, and the shape of the 
vortex near the wall is the same as in the bulk. Thus, an initial 
quarter ring vortex with a radius small compared to the cylin- 
der radius is closer to the non-oscillating configuration than a 
quarter ring vortex with a radius comparable to the cylinder 
radius. 



VI. EFFECT OF TILTING 

Tilting the vessel has the advantage of breaking the cylin- 
drical symmetry. A small tilt may be useful in calculations 
also for detecting effects that occur due to this asymmetry, 
unavoidably present in experimental situations. 

The most noticeable effect of tilting is the oscillation of vl z 
about a constant value (close to aQR). This is due to a sinu- 
soidal component of the boundary field 16 . 

A somewhat surprising effect was also found. At a large 
enough tilt angle the azimuthal velocity of the vortex end may 
approach to zero (in moving coordinates). Thus, in the asymp- 
totic limit, the vortex end becomes locked at some azimuthal 
angle, while a constant axial velocity component remains. A 
detailed study of this phenomenon has not been done yet, but 
this sort of behavior can be seen in the last two pictures of 
Fig. 10, corresponding to tilting angles 50° and 60°. 



scale vortex tangle calculations to do a very simple job, there 
still is some potentially useful information to be gained from 
this endeavor. 

Vortex motion was found to be quite stable in both the 
untilted and moderately tilted (< 30°) cases. This lim- 
its the causes for the experimentally observed single vortex 
instability 2 ' 11 ' 14 to areas not covered in this study, such as sur- 
face effects and dissipative bulk effects. 

The asymptotic line velocities aQR and (a' — 1)QR, valid 
in the flR 2 / k ^> 1 limit for the untilted case, could be gen- 
eralized by equations of type (11) with some yet unexplained 
numerical constants. Providing an explanation for these con- 
stants might necessitate analyzing, either analytically or nu- 
merically, the shape of the vortex, or including varying also 
the parameter R/a. 

This study also highlights the role of scaling properties in 
superfluids and other analogous systems. A deeper analysis 
to connect these formulas to the velocity of the vortex front 
should be carried out. The azimuthally locked motion in the 
tilted case was an unexpected result. 
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